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Kinetic equation with exact charge conservation 
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We formulate the kinetic master equation describing the production of charged particles which are 
created or destroyed only in pairs due to the conservation of their Abelian charge. Our equation 
applies to arbitrary particle multiplicities and reproduces the equilibrium results for both canoni- 
cal (rare particles) and grand canonical (abundant particles) systems. For canonical systems, the 
equilibrium multiplicity is much lower and the relaxation time is much shorter than the naive extrap- 
olation from the grand canonical ensemble results. Implications for particle chemical equilibration 
in heavy-ion collisions are discussed. 
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I. INTRODUCTION 

Relativistic statistical thermodynamics has long been 
used as a tool to describe particle production in heavy- 
ion and in high-energy particle collisions flpf|- Recent 
analyses have shown that the statistical models can in- 
deed give a satisfactory description of the multiplicities 
of most hadrons measured in A-A collisions at AGS and 
SPS energies However, the dynamics of particle 

equilibration, and in particular chemical equilibration, is 
still not well understood In this Letter, we shall 

address the problem of chemical equilibration within the 
statistical kinetic approach. 

Within this approach, particle production is commonly 
described using the grand canonical ensemble, where 
event-averaged multiplicities are controlled by chemical 
potentials. In this description the net value of a given 
17(1) charge (e.g., electric, baryon, strangeness, charm, 
etc.) fluctuates. These fluctuations can be neglected 
only if the particles carrying the charge in question are 
abundant. In this case the charge will be conserved on 
the average. In the opposite limit of rare particle pro- 
duction, conservation laws must be implemented locally 
on an event-by-event basis |fl,p|Js|Jlot| ; i.e., a canonical en- 
semble must be used. 

The local conservation of quantum numbers in the 
canonical approach severely reduces the phase space 
available for particle production |l|J^,^|,p|-^2| . Recently, 
it has been shown that the canonical statistical model 
provides a good description of particle yields measured 
in low-energy heavy-ion [p"3[ and high-energy hadron- 
nucleus, hadron-hadron and e + e~ reactions [pt pd|JT5[| . 

In view of these results it is of particular importance 
to formulate a kinetic theory for the time evolution of 



particle production in order to investigate the approach 
to the canonical chemical equilibrium. Obviously, if par- 
ticles are abundantly produced, the equilibrium result for 
the particle multiplicity should coincide with the grand 
canonical one. For rare processes, however, the parti- 
cle production is strongly correlated and the canonical 
equilibrium result is expected. The kinetic formulation 
for the production of strongly correlated particles was 
studied in the literature However, no complete 

solution has been obtained. 

In this letter we consider the time evolution of the mul- 
tiplicity of particles that carry the charges corresponding 
to an U(l) internal symmetry. We formulate a kinetic 
master equation valid for arbitrary particle multiplicity. 
It reproduces the canonical equilibrium solution for rare 
particle production, and reduces to the standard grand 
canonical rate equation for abundant particle production. 



II. RATE EQUATION 



In the standard formulation [ |16|Jlq -|2fj[], the rate equa- 
tion for a binary process a\a,2 — > 01&2 with a 7^ b is 
described by the following population equation: 



(1) 



where G = {(Jqv) and L = (<tl,v) give the momentum- 
averaged cross sections for the gain process 040,2 — > &162 
and the loss process 6162 — * 0402, respectively. AT*, repre- 
sents the total number of particles k, and V is the proper 
volume. Among such processes, a typical example is the 
kaon production/annihilation via 7r + 7r~ <-> K + K~ . The 
above rate equation, however, cannot be applied to the 
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situation where particle production is rare and is strongly 
correlated by exact charge conservation. 

To account for the correlation between the produc- 
tion/annihilation of particles b\ and b 2 , let us define Pij 
as the probability to find a number i of particle b\ and 
a number j of particle b 2 in an event. We also denote 
as Pi the probability to find a number i of particle b in 
an event. The average number of particles b per event is 
then defined as: 



(N b 



IP,, 



(2) 



We can now write the following general rate equation 
for the average particle multiplicity: 



(3) 



Further, let particles b\ and b 2 carry opposite units of 
a charge, corresponding to an U(l) internal symmetry 
(strangeness in the case of kaons). Then the U(l) charge 
neutrality of the system gives TV = N bl = N b2 . We have 
then, 

p — p. x. . 

£ ijP itj = £ ^ 2 J R i = (TV 2 ) = (TV) 2 + (SN 2 ), (4) 

i,j » 

where (SN 2 ) represents the event- by- event fluctuation of 
the number of 6162 pairs. Note that we always consider a± 
and a 2 particles abundant (such as, e.g., pions) so that 
we can neglect the event- by- event fluctuations of their 
multiplicity and the change of their number due to the 
considered processes. 

Following Eqs.(|3|) and (Q), the general rate equation 
for the average number of b\b 2 pairs can be written as 



(5) 



For abundant production of b\b 2 pairs, where (TV) 1, 
(TV 2 ) « (TV) 2 , (6) 
and Eq.(^) obviously reduces to the standard form, i.e., 

(7) 



d(N) G, i; , L, ,, 



However, for rare production of b±b 2 pairs, where (TV) <C 
1, the rate equations ((TJ) and (0) are no longer valid. We 
have instead 



(TV 2 )«(TV), 
which reduces Eq.(^J) to the following form: 

Thus, in the limit (TV) -c 1, the absorption term depends 
on the pair number only linearly, instead of quadratically 
for the limit (TV) > 1. 



(8) 



(9) 



III. EQUILIBRIUM MULTIPLICITY AND 
RELAXATION TIME 

To illustrate the differences in the time evolution of 
particle abundance, we consider the two limiting cases, 
(TV) ^> 1 and (TV) <C 1, and present their equilibrium 
values and relaxation times for the production of 6162 
pairs. As an example, we consider a system at fixed 
temperature and volume and with no initial 6162 pairs, 
i.e., (TV)(r = 0) = 0. 

In the limit when (TV) 3> 1, the standard Eq.(Q) is valid 
and has the following well-known solution: 



(N) GC (t)= TV c G q c tanh (r/r GC ) 



(10) 



where the equilibrium value for the number of &1&2 pairs 
TV^ C and the relaxation time constant t^ c are given by 



GC 




V 
T7e 



(11) 



respectively, with e = G(N ai )(N a2 ) /L. 

In the particular case where particle momentum dis- 
tributions are thermal, the gain (G) and loss (L) terms 
just represent the thermal averages of the production and 
absorption cross sections with 



G _ d hl al i K2(a bl )d b2 al 2 K 2 (a b2 ) 
L d ai a 2 . i K 2 (a ai )d 



(12) 



where dfc's denote the degeneracy factors, and ctk = 
rrik/T. The equilibrium value for the number of 6162 
pairs in Eq. (|l l[) now reads as 



TV GC 



^VT 3 a 2 bl K 2 (a bl ] 



(13) 



Thus it is described by the Grand Canonical (GC) result 
with vanishing chemical potential due to our requirement 
of the (average) U(l) charge neutrality of the system. 

In the opposite limit where (TV) -C 1, the time evo- 
lution is described by Eq.(^J), which has the following 
solution: 



(N) c (T) = NC q (l-e-r^y (14) 
with the equilibrium value and relaxation time given by 

(15) 



TV° = f 



V 
L' 



With a thermal momentum distribution the equilib- 
rium value of 6162 pair multiplicity has the following 
form: 



TV C = 

1 *cq 



^VT*alK 2 (a bl ) 



^VT*a 2 b2 K 2 (a b2 



(16) 



This equation demonstrates the locality of the U(l) 
charge conservation. With each particle 61 , a particle b 2 
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with the opposite charge is produced in the same event 
in order to conserve charge locally. This is the result 
expected from the Canonical (C) formulation of conser- 
vation laws 

We note that Eq. (jilt) is just the leading term in the 
expansion of the canonical result for the multiplicity of 
particles carrying the U(l) charges. The general expres- 
sion is known to have the following form |jl(|fi"i"|| : 



Nt. 



N, 



cc * I (2Ngcy 



GO 



(17) 



where A^ c is given by Eq.(|T|) and /, 's are the modified 
Bessel functions. 

Comparing Eq. (|TT|) and Eq.([l5|), we first find that, for 
(N) <C 1, the equilibrium multiplicity in the canonical 
formulation is much lower than what is expected from 
the grand canonical result, 



(iV G q c ) 2 « TV-gc. 



(18) 



This shows the importance of the canonical description 
of charge conservation when the multiplicity of charged 
particles is small. We also note that the volume depen- 
dence in the two cases is different. The particle density 
in the GC (abundant) limit is independent of V, whereas 
in the opposite canonical (rare) limit the density scales 
linearly with V . 

Secondly, the relaxation time for a canonical system is 
far shorter than what is expected from the grand canon- 
ical result, 



'0 



cq 



GC 



(19) 



due to small number of particles (iV^ c <C 1). For ex- 
ample, the total number of produced kaons in Au + Au 
collisions at 1 GeV/A is of the order of 0.02. Thus the 
canonical relaxation time is a factor of 50 shorter than 
what is expected from the grand canonical formulation. 

We note from Eq.(||) that these two limits are essen- 
tially determined by the size of (SN 2 ), the event-by- 
event fluctuation of the number of &1&2 pairs. The grand 
canonical results correspond to small fluctuations, i.e., 
(8N 2 ) I (N) 2 << 1, while the canonical description is nec- 
essary in the opposite limit. 



IV. MASTER EQUATION 

In this section, we formulate the general evolution 
equation which is valid for an arbitrary value of (N). 
It is a master equation for P n (r) , the probability of find- 
ing n pairs of &1&2 at time r. This probability increases 
with time due to transitions from n — 1 and n + 1 states 
to the n state, while it also decreases due to transitions 
from the n state to n— 1 and n + 1 states. The transition 
probability n — ► n+1 per unit time due to pair creation is 



G{N ai ){Na 2 ) /V and the transition probability n — > n— 1 
due to pair annihilation is n 2 L/V. Therefore, the master 
equation set has the form: 



Clin Cx , , , 



,) [Pn-l 
(n + 1) 2 P„+ 



Pn] 



(20) 



where n = 0, 1, 2, 3, • • •. Multiplying the above equation 
by n and summing over n, one recovers Eq.(^J), the gen- 
eral rate equation for the time evolution of the average 
number of 6162 pairs. However, the master equation j20| ) 
contains much more information than the rate equation 
(||). It contains enough information to solve for the evo- 
lution of (N) (and all moments of TV) for arbitrary (TV). 
For example, one can also obtain the time evolution of 
particle fluctuations, (SN 2 ), which are of physical impor- 
tance for rarely produced particles. 

We can convert the iterative equations ( p0| ) for P n 's 
into a partial differential equation for the generating 
function 



00 

3 (:r,T) = ]>>"P n (T). 

n=0 



(21) 



Multiplying Eq.(|20|) by x n and summing over n, we find 
dg(x,r) L 



:(l-x) (xg" + g'-eg), 



(22) 



dr V 

where g' = dg/dx. Note that the g(l, r) does not change 
with time, which is equivalent to the conservation of total 
probability evident in Eq.(|2fj|). 

The equilibrium solution, g cq (x), thus, obeys the fol- 
lowing equation: 

X 9eq + 9cq ~ e 3eq = 0. (23) 

By variable substitution (x = y 2 /(4e)) this equation can 
be reduced to the Bessel equation. The solution that is 
regular at x = (since g(0) = P < 1) is given by 
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(24) 



where the normalization is fixed byg(l) = P n = 1. 

The equilibrium probability distribution P n can now 
be found from Eqs.(ElLE4h as 



Pn 



7 (2Vi) (n!) 2 ' 



(25) 



and the value of the average number of &162 pairs per 
event at equilibrium is given by 



(iV) eq = g\l) = V~e 



MVg) 

/o(2v^)' 



(26) 



which obviously coincides with the expected result for 
particle multiplicity in the canonical ensemble given by 
Eq.(|l7|) and also reduces to the results for the two limit- 
ing cases in Sec. HI. 
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V. CONCLUSIONS 

We have formulated the kinetic master equation for 
strongly correlated production of particles, where the cor- 
relation is due to the local charge conservation required 
by an U(l) internal symmetry. Our general rate equation 
is valid for arbitrary value of (N) , thus it reduces to the 
grand canonical results for large (N) and to the canonical 
results for small (N). Therefore, our equation provides 
a generalization of the standard rate equation beyond 
the grand canonical limit. We have shown that for rare 
particle production the equilibrium multiplicity is much 
smaller and the relaxation time is much shorter than ex- 
pected from the standard rate equation. For abundant 
particle production, where the standard rate equation ap- 
plies to first order, one can use the general rate equations 
to study finite-size corrections to the grand canonical re- 
sults. 

Our results could be of importance in the description 
and understanding of equilibration phenomena and equi- 
librium properties of partonic or hadronic matter cre- 
ated in heavy-ion collisions. For example, it could pro- 
vide insights into the equilibration time of strange parti- 
cles produced at SIS, or open charm and charmonium 
productions at SPS and higher energies. It may also 
be meaningful for transport model studies of rare par- 
ticle production, especially the perturbative procedure 
for rare processes pH , where the local charge conserva- 
tion requires that the production probability should be 
assigned to the pair instead to each particle separately. 
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